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ABSTRACT 

The possibility to unambiguously determine the equation-of-state of the cosmic dark 
energy with existing and future supernovae data is investigated. We consider four 
evolution laws for this equation-of-state corresponding to four quintessential models, 
i.e. i) a cosmological constant, ii) a general barotropic fluid, iii) a perfect fluid with 
a linear equation-of-state and iv) a more physical model based on a pseudo-Nambu- 
Goldstone boson field. We explicitly show the degeneracies present not only within each 
model but also between the different models : they are caused by the multi-integral 
relation between the equation-of-state of dark energy and the luminosity distance. 
Present supernova observations are analysed using a standard method and the 
minimal values obtained for each model are compared. We confirm the difficulty 
to discriminate between these models using present SNela data only. By means of 
simulations, we then show that future SNAP observations will not remove all the 
degeneracies. For example, wrong estimations of Q,m with a good value of Xmin could 
be found if the right cosmological model is not used to fit the data. We finally give some 
probabilities to obtain unambiguous results, free from degeneracies. In particular, the 
probability to confuse a cosmological constant with a true barotropic fluid with an 
equation-of-state different from —1 is shown to be 95 % at a 2 cr level. 

Key words: cosmological parameters - dark matter - supernovae : general 



1 INTRODUCTION 

Present supernovae data strongly support cosmological 
models containing a perfect fluid with a negative pressure 
(Riess et al. 1998, Perlmutter et al. 1999). The oldest and 
most studied candidate for this fluid is the cosmological con- 
stant, which acts like a perfect fluid whose equation-of-state 
(hereafter, EOS) is w = p / p = —1, where p is the fluid pres- 
sure and p its density. But the vacuum energy density associ- 
ated with the cosmological constant is 60-120 orders of mag- 
nitude smaller than its natural value derived from quantum 
field theories. This discrepancy is known as the cosmologi- 
cal constant problem (Abbott 1988, Weinberg 1989, Carroll 
et al. 1992, Sahni & Starobinsky 2000) and has led theorists 
to find alternative dark energy candidates with a present 
negative pressure. The simplest models are based on a gen- 
eralisation of the cosmological constant, e.g. a barotropic 
fiuid {w{z) — Wo = constant < 1 ; Gonzalez-Diaz 2000, Di 
Pietro & Demaret 2001) or a homogeneous fiuid for which 
ui is a linear function of z {w{z) = wo + wi z with wo and 
wi constant; Goliath et al. 2001, Maor et al. 2001, 2002). 

These models are interesting because they are described 
by simple field equations. However, they suffer from a lack 
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of physical justification. Other models with a more gen- 
eral EOS and stronger physical interpretation have been 
proposed (Peebles & Ratra 1988, Ratra & Peebles 1988, 
Wetterich 1988, Ferreira & Joyce 1998, Steinhardt et al. 
1999). The most popular ones are based on a dynamical 
quintessential component represented by a minimally cou- 
pled scalar field evolving in a potential. Mathematically, 
such a quintessence fluid can be described by an EOS w 
function of the redshift 2 (— 1 ^ w{z) ^ 1). Among these 
quintessence models, we shall consider the one which as- 
sumes the existence of an ultra-light pseudo-Nambu Gold- 
stone boson (PNGB) field relaxing to its vacuum state (Frie- 
man & Waga 1998, Waga & Frieman 2000, Ng & Wiffshire 
2001a). From the quantum viewpoint, the PNGB models are 
the simplest way to introduce a naturally ultra-light scalar 
field able to reproduce the cosmological observations (Frie- 
man et al. 1992, 1995). Moreover PNGB models provide an 
interesting theoretical framework to any spontaneous sym- 
metry breaking which could justify the neutrino mass found 
in the Mikheyev-Smirnov-Wolfenstein solution to the solar 
neutrino problem (Wolfenstein 1979, Mikheyev & Smirnov 
1985). 

In summary, many theoretical models describing a cos- 
mological fiuid with negative pressure have been proposed : 
some are designed to be mathematically simple while others 
rely on a more physical justification. The first aim of this pa- 
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per is to put in evidence the strong degeneracies existing be- 
tween the luminosity distance predicted by four quintessence 
models : a cosmological constant, a general barotropic fluid, 
a fluid with a linear equation-of-state and a more physi- 
cal model based on a pseudo-Nambu-Goldstone boson field. 
Second, we shall determine the constraints on those models 
coming from the present observations of type la supernovae 
(hereafter, SNela) and confirm the need for more data to 
discriminate between them. This may be fulfilled by the pro- 
posed SNAP satellite (SuperNova/Acceleration Probe ; see 
SNAP URL). The main objective of this instrument is to 
detect a very large number of supernovae up to a redshift 
of 1.7, in order to yield a more precise determination of the 
cosmological parameters and therefore to provide informa- 
tion on the nature of dark energy. The third objective of 
this paper is to analyse, by means of simulated data, how 
well future SNAP observations alone will be able to break 
the degeneracies that we put in evidence. 

Several authors have already studied the feasibility of 
SNAP to determine the properties of the dark energy. De- 
pending on the method used for reconstructing the cosmo- 
logical model, their conclusions are quite different : some au- 
thors are optimistic regarding the possible determination of 
the EOS of the dark energy using SNela (Huterer & Turner 
1999, Nakaruma & Chiba 1999, Saini et al 2000, Chiba & 
Nakaruma 2000, Weller & Albrecht 2001, 2002) while oth- 
ers are more cautious (Barger & Marfatia 2001, Astier 2001, 
Maor et al. 2001, 2002, Gerke & Efstathiou 2002). Most of 
these discrepancies can be traced back to the differences in 
initial assumptions, prior knowledge, ... used for the recon- 
struction of the cosmological model (Goliath 2001). Usually, 
optimistic conclusions result from strong assumptions, such 
as an accurate knowledge of i},n. 

The structure of the paper is as follows. In Section 2, 
we present the field equations describing the four cosmolog- 
ical models we have considered. In Section 3, we explicitly 
display the degeneracies between the luminosity distances 
predicted by these four models. The constraints brought by 
present SNela data on the parameters of each model are 
presented in Section 4. In Section 5, we explore the ability 
of SNAP to break the degeneracies and so to discriminate 
among these models. Finally, we summarize our work and 
give some conclusions in Section 6. 



2 QUINTESSENCE MODELS : THEORETICAL 
FRAMEWORK 

The field equations of a Friedmann-Lemaltre-Robertson- 
Walker spacetime filled with two non-interacting fiuids, a 
pressureless fiuid (matter - flm) and a spatially homoge- 
neous scalar field </> in a potential V{(f>) (quintessence - Qx), 
can be written as 



Hi 



f^m (1 + zf + f^fc (1 + zf + Qx g{z) 



Px{z) = ~3[pxiz) +px(z)]H{z) 
= ^3[1 + w{z)]p4z)H{z), 



(1) 



(2) 



where the prime denotes the time derivative, H{z), the 
Hubble parameter, Qk, the curvature parameter {Qk = 
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with K = 8 TT G. 

When the scalar field acts like a fiuid with a linear EOS 
in z, i.e. 'w{z) = wq -\- wi z, equation JTJ becomes g[z) = 
g3»i z _^^^3{i-^»o-»i) and the Friedmann equation given 
by Q can be written as 

H^{z) = Hl{{l + zf[^mz+l] 

-f nx(l + zf [e^*"!" (1 + 2)l + 3(»0-»l) _ I ^ (g) 

where the relation between the density parameters has been 
used, i.e. fit = 1 — Q,m — ^x- 

The luminosity distance dL{z) to an object at redshift 
z is such that 



Dt{z) 



dLiz)Ho 




(9) 
(10) 



with S{x) = sin(a::), x, sinh(2;) for respectively a spatially 
closed, flat and open universe. Note that di is the luminosity 
distance whereas is a dimensionless quantity. For in 
units of Mpc, the apparent magnitude can be expressed in 
terms of D^j^ as follows 



niB = 5 log(_Di 



Me 



(11) 



where Mb is the magnitude zero point in the B band defined 
by 



Mb = Mb ~ 5 log[Ho/c]+ 25 



(12) 



with Mb , the absolute magnitude of the SNela. 

Hereafter, we focus on three types of models containing 
a scalar fluid with a linear EOS (cf. also models A, B and C 
in Table [U : 

- Model A : a scalar field still acting as a cosmological 
constant with no constraint on the geometry of the universe, 
i. e. a FLRW model, 

- Model B : a scalar fluid with a constant EOS in a spa- 
tially fiat universe, 

- Model C : a scalar fiuid with a hnear EOS in a spatially 
fiat universe. 

These models are mathematically the simplest ones but 
they have no satisfactory physical interpretation. So we shall 
also consider a fourth quintessence model, based on phys- 
ical arguments. This model is constructed with an ultra- 
light pseudo-Nambu-Goldstone boson (PNGB) field (Frie- 
man & Waga 1998, Waga & Frieman 2000, Ng & Wiltshire 
2001a, 2001b). In quantum field theory, a Nambu-Goldstone 
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Name 


Constraints 




Parameters 


Model A 


wq = — 1 and wi = 






Model B 


Qm + f^x = 1 and wi 


= 


^m, WO 


Model C 


Qrn + = 1 and S7,„ 


= 0.3 


Wo, Wl 


Model D 


Qrn + ^x = 1 and 




M, f 




V'(</.) = M4 {cos[<A//] 


+ 1) 





Table 1 . Characteristics of the four theoretical models considered 
in this paper. 



boson field is generated wlien a global symmetry is spon- 
taneously broken. Moreover, if this symmetry is explicitly 
broken, the Nambu-Goldstone boson acquires a mass and 
becomes a pseudo-Nambu-Goldstone boson. A quintessence 
model based on a PNGB field can be described by the field 
equations (0-||7|l with the following periodic potential (Frie- 
man et al. 1992, 1995) : 



=M* (cos [^/f] + l) 



(13) 



where M and / are the two mass scales, i.e. / is the 
spontaneous- and M is the explicit symmetry breaking scale. 
For such a model, the field equations Q-Q have no analyt- 
ical solution. Therefore, the behavior of the functions g(z) 
and H{z) can only be derived from a numerical integration. 
We take as initial conditions (f>[ — and (f>i — 10^^ mpi 
where mpi is the reduced Planck mass (see justification be- 
low). Assuming again a spatially fiat universe, we are left 
with a model depending on only two parameters, M and /. 
This model is referred to as " Model D " (cf. also Table 
and its most important properties are summarized in Fig- 
ures Q and |5] Because of the Hubble damping, the PNGB 
field is initially frozen near the top of its potential ( ~ 
and -C /) and therefore acts as a cosmological constant, 
so that w{z ^ 1) « —1. When the universe temperature 
becomes less than the PNGB mass, i.e. — / f , the 
field starts to slowly roll down the potential. The minimum 
of the potential is reached asymptotically, after a large num- 
ber of coherent oscillations about it. This second phase cor- 
responds to a scalar component which behaves on average 
as nonrelativistic matter. It means that model D looks ini- 
tially like model A with ilk = and asymptotically like 
an Einstein-de Sitter universe. The larger the quintessence 
mass, the sooner the asymptotic state is reached. If the 
scalar field evolution starts before its energy density becomes 
dominant, the PNGB field will never be relevant for the uni- 
verse dynamics (see Figure 0\.). However, if the dynamical 
phase of the scalar field begins when the vacuum energy as- 
sociated with the PNGB field is already dominant, then the 
scalar energy density may dominate today the cosmic den- 
sity of the universe and provide a non-negligible dark energy 
component (see Figure 03). 

Figure |21 shows examples of curves of equal Qm in an 
M — f diagram. The shaded area represents the parame- 
ter regions in which the expansion is accelerating today. In 
these regions, the kinetic energy of the PNGB field is small 
compared to its potential. This occurs either when the scalar 
field is still frozen (the lowest and largest shaded band) or 
when its first derivative changes sign, i.e. when <j)' ^ (the 
thin shaded bands). Regions outside these shaded areas cor- 
respond to models where the scalar field is presently near the 





Figure 1. Behavior of (thick) and pm (thin), in units of the 
present critical density, when the quintessence component is a 
PNGB field. In the small boxes, the corresponding behavior of 
ui0 is given. Negative redshifts are associated with the future of 
= 1.4 X lO-^hj/^ eV and / = 3 X 10^^ 
4 X lO-'^h}/'^ eV and f = 5 X lO^'^ GeV. 
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Figure 2. M — f diagram for the PNGB model. The shaded areas 
represent the parameter regions where the universe is presently 
accelerating whereas the solid lines are the curves of equal f2„i 
with Qm = 0.8, 0.6, 0.4, 0.2, from left to right. 



minimum of its potential (for more details on PNGB mod- 
els, see Frieman et al. 1992, 1995, Frieman & Waga 1998, 
Waga & Frieman 2000, Ng & Wihshire 2001a). 



3 DEGENERACIES BETWEEN THE 
QUINTESSENCE MODELS 

Since the SNela apparent magnitudes are sensitive to the 
luminosity distance and since the latter is related to the 
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Figure 3. Curves of equal luminosity distance in the framework of the four cosmological models considered in this paper, (a) Model A : the 
dashed curves represent the set of values for {Qm, ^a) such that Di^(zi; flm, ^a) = Di^{zi\ 0, 0), with Zi = 0.05, 0.4, 0.8, 1.0, 1.2, 1.4. Sim- 
ilarly, the solid curves show the set of values for (n,„, VLa) such that DL{zi;Q,m, S^a) = (zi; 0, — 1), with Zi = 0.05, 0.4, 0.8, 1.0, 1.2, 1.4. 
(b) Model B : the curves represent the set of values for (Hm. mn) such that Djjzi; Vlm.wn) = Di^(zi; 0.7, —0.5) (top) and ^^(zi; Cm, two) = 
Dl (zi; 0.3, —0.7) (bottom), with Zi = 0.05, 0.4, 0.8, 1, 1.2, 1.4. (c) Model C : the curves represent the set of values for {wq,wi) such that 
DL{zi; wo, wi) = D£(zj;-0.3,0.2) (top) and Di(zi; too, = Di(zi; -0.7, -0.5) (bottom), with zt = 0.05,0.4,0.8,1,1.2 and 1.4. (d) 
Model D : curves show the set of values of (/, Af) such that Di(zi;/, M) = Di(zi;0.5,3) with Zi = 0.4 (thin), 0.8 (thick) and 1.2 
(dashed). 



model parameters through a multiple integral (see eq. IIOII 
with (0 and 0), it is quite natural to expect some de- 
generacies between the different models. This is even true 
within a given cosmoiogicai modef. Indeed, Figure |21 shows 
curves of equal luminosity distance for different redshifts, in 
the parameter spaces of the four models. For each model, 
the curve shape varies with the given redshift. The curves 
intersect in one point in models A, B and C and this lifts 
the degeneracy. But in model D, the contours could remain 
close to each other, making more difficult a clear estimation 
of the cosmological parameters. 

The situation is even worse when comparing the dif- 
ferent models. Indeed, Figure |1] shows expected magnitude 
differences between different models, versus redshift. Spe- 
cific choices of parameters of either model A, B, C or D can 
yield magnitude differences less than 0.04 mag all the way 
to z = 2. These are really small compared to the 0.15 mag 
intrinsic dispersion of the SNela maximum luminosity. In 
other words, as already mentionned by Maor et al. (2001, 
2002), different function 'w{z) can approximatively lead to 
the same dL{z) relation. A large number of SNela obser- 
vations with very good photometry are needed in order to 



break these degeneracies. At this point, data at large red- 
shifts {z > 2) could be more than useful. 

Figure |5] illustrates that estimates of the cosmological 
parameters can simply be erroneous if a wrong cosmological 
model is chosen for the analysis. Here are the different steps 
followed to draw this plot : 

★ We chose some values of the parameters describing the 
model C, i.e. (wo, wi) = (—0.7, —0.5). 

★ The luminosity distance predicted by model C at z = 
0.05 were numerically computed. 

★ We then searched the parameters {flm, wo) for which 
model B predicts the same luminosity distance at z = 0.05. 
We plotted the curve of equal luminosity distance obtained 
in the parameter space of model B. 

★ We repeated the two previous steps for different val- 
ues of z, i.e. z = 0.4, 0.8, 1, 1.2, 1.4, 1.7, 3, and plotted the 6 
corresponding curves. 

★ The four previous operations were repeated with a sec- 
ond set of parameters for the model C, i.e. {wo,wi) = 
(-0.3,0.2). 

Keeping in mind that models C are defined with f2m ~ 
0.3, Figure 5 clearly shows that the value of i},n would not 
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Figure 4. Comparison between the magnitudes pre(iicted by var- 
ious models A and B (a), models A and C (h), models B and C 
(c) and models C and D (d). 



be properly recovered by fitting model B on the data. In- 
deed, the first set of curves ( correspondong to (wojWi) = 
(—0.7, —0.5) in model C) all cross over a very small area 
characterized by (fim, wo) ~ (0.15, —0.6). Similarly, the sec- 
ond set of curves (corresponding to {wq,w\) = (—0.3,0.2) 
in model C) would lead to an estimate of (f2m,'ii'o) ~ 
(0.65, —0.6) if the data analysis was carried out with model 
B. These results mean that the same data can lead to differ- 
ent estimates of f2m and wa, depending on the cosmological 
model used for the data analysis. In other words, the role 
of the derivative of the dark energy EOS on the luminosity 
distance can be played by a matter density parameter in 
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Figure 5. Curves in the parameter space of model B 
such that Dj^{zi\Q.m,WQ) built with this model is equal to 
Di^{zi\WQ = — 0.7,toi = —0.5) (left) and to D]:^{zi\wo = 
— 0.3,uii = 0.2) (right) derived from model C , for Zi = 
0.05, 0.4, 0.8, 1, 1.2, 1.4, 1.7, 3. The presence of degeneracies is re- 
vealed by the fact that the different curves cross each other over 
a very small area in the parameter space. The curves correspond- 
ing to 2 = 3 (dotted) indicate that the degeneracy could break at 
very large redshifts. 



models with a constant EOS : a model with a fixed Q.m con- 
taining a fluid described by a negative (positive) w\ can be 
mistaken with a model described by a smaller (larger) Q.m 
containing a barotropic fluid (see also Figure [T^ . Note that 
the dotted curve corresponding to z = 2> shows degeneracy 
breaking at very large redshifts (see also Figure |4j. 



4 QUINTESSENCE MODELS AND PRESENT 
SNEIA DATA 

We used the published data of Perlmutter et al. (1999; P99 
hereafter), which consist of a sample containing 60 SNela, 
to determine the constraints that we have today on the pa- 
rameters describing the four models introduced in Section 2. 
As it has been noted in P99, four of these supernovae are 
"outliers" . So we excluded them and considered a sample of 
only 56 SNela^ . We analysed these data with a standard 
method, i.e. by minimizing the value of defined by 



56 

where m\ 



(14) 



(resp. rrig ") are the effective apparent mag- 
nitudes (i. e. the observed magnitudes corrected by the 
lightcurve width-luminosity relation, the Galactic extinc- 
tion, the K-correction and the cross-filter calibration) pre- 
dicted by the model (resp. given by the data), ai are the 
uncertainties on the data and 9i et 62 are the two parame- 
ters of the theoretical model (cf. Table 0. 



^ Note that using the sample of 60 SNela does not significantly 
change the results: the contour positions are approximately the 
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Figure 6. 68.3% (dark) and 95.4% (light) confidence contours 
obtained by analysing P99 data with each model. Solid curves 
separate universes which are presently decelerating from those 
which are presently accelerating. The four plots on the left (resp. 
on the right) are obtained with Mg = 24 (resp. with no prior 
constraints on Mb)- 



Concerning Mb, we adopted two different approaches. 
We first assumed that Mb was exactly known and equal to 



24,^ all uncertainties on this parameter being neglected. In 
the second approach, we assumed no prior constraint on Mb , 
which is just an unknown constant with a value between — oo 
and +00. In this case, we integrated the probabilities on Mb 



and therefore worked with a x defined by 



X = -2 In 



dMB = 



where 

^ _ [mB.i -5 1og[I)^]j 

1—1 



B = 



C 



E 

j=i 

56 



ml^f - 5 log[Di.'^ 



(16) 
(17) 
(18) 



Note that a uniform distribution of the probabilities on a 
relatively small interval, i.e. Mb = 24 ± 0.2, gives the same 
results as a uniform distribution between — oo and +oo. 

We analysed the sample of 56 SNela given in P99 us- 
ing as framework the four theoretical models discussed in 
Section 2 and for two opposite assumptions on Mb- The 
probability contours are shown in Figure|S| Table|5|lists the 
best fit values of the parameters together with the Xmin a^nd 
Xmin/DOF (DOF: Degree of freedom). 

The results obtained using model A are in agreement 
with those found in the literature (P99). Small differences 
may be related to the fact that we used data already cor- 
rected for the stretch-factor whereas P99 integrate on this 
parameter. However, these differences are not relevant in 
view of our objectives. Results concerning model B show 
that present data agree with a baroropic fluid characterized 
by wq < —0.4 at the 2 a level. An important remark has to 
be made concerning models B and C : FigureElindicates that 
the data favour unphysical regions where wo < — 1 (see also 
Caldwell 2002). Note that if we restrict ourselves to mod- 
els with Wo ^ —1, the best fit positions, characterized by 
(ilm. Wo) = (0.3, —1) in model B and by {wo,wi) = (—1, 0.6) 
in model C, are associated to Xmin values only marginally 
different from those obtained when wo < —1 (Axmin < 0.1). 
The fact that the true minimum is located in a unphysical 
region is something that has to be strongly considered in the 
future because, if confirmed, it can rule out these models. 

The constraints for model D select two sorts of regions 
(Figure il; see also Ng & Wihshire 2001b) : one where the 
scalar field is still frozen (lower horizontal band) and another 
where its evolution has just begun (upper oblique band). 
Note that these contours are not closed : the larger the /, 
the more the two branches are separated. A small part of 
the 2 cr contours also appears on the second plot of model D, 
selecting a region with M ~ 6 X 10"^ hl^^ eV. If we adopt 
the point-of-view that the PNGB field of model D is related 
to the resolution of the solar neutrino problem, then the 
value of M gives an estimate of the neutrino mass (Frieman 



^ This value is the one obtained when Mb = —19.5 mag and 
same and the corresponding P^r degree of freedom is slightly Hq = 58.5 km/s/Mpc are introduced in relation I12i (Parodi et 

larger. al. 2000). 
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= 0.6 


Qa = 0.9 


B 


62.2 


1.15 




= 0.41 


wo = -1.2 


C 


62.2 


1.15 


WQ -- 


= -1.1 


«)i = 0.9 
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62.2 


1.15 


M 


= 2.3 


/ = 0.5 




Model 


■ 










A 


66.8 


1.26 




= 0.8 


Ca = 1.3 


B 


67.4 


1.27 


S^m 


= 0.46 


wo = —1.6 


C 


67.3 


1.27 


wo = 


= -1.4 


= 2.3 


D 


67.4 


1.27 


M 


= 2.3 


/ = 0.5 



Table 2. Minimum values of a-nd obtained when the 
analysis of the existing data is carried out on each of the four 
quintessence models. The upper table assumes Mg = 24 whereas 
the lower one assumes no prior constraints on Mb ■ 



et al. 1992, 1995). The last two plots of Figure El show that 
at least one neutrino species has a mass of: ~ 2 x 10~^ ftj''^ 
eV, if the PNGB field is still frozen, ~ 3 x 10"^ hl'^ eV, if 
the field is just starting to evolve or ~ 6 x 10~^ /ij''^ eV (2 a 
region on the second plot), if the field has already passed 
one time through the potential minimum. These estimates 
are indeed similar to the ones found by SuperKamiokande 
and SNO (Sudbury Neutrino Observatory) for the neutrino 
mass squared differences, i.e. 5m\-j — (1.6 — 3.9) x 10~^ 
eV^ (Fukuda et al. 1998) and (Jmfa = (0.3 - 3.5) x IQ-" 
eV'^ (Choubey et al. 2002 and references therein) . This may 
be only a numerical coincidence... or it may comfort us in 
the resolution of the neutrino problem through a Goldstone 
mechanism. 

Note that contours where M --^ 3 x 10""^ eV or 
are associated to small values of Q.rn (<?C 0.1 ; 
see Figurel^l which are inconsistent with observations (Pee- 
bles & Ratra 2002) . The value of / is only poorly constrained 
by the present data. This is because / becomes an impor- 
tant parameter for the cosmological model only when the 
scalar field evolves in its potential, not when it is frozen or 
just starting to move. 

As it can be seen on Figure |S] in some cases, present 
data can also favour a decelerating universe, even at 1 a 
level. An important remark has to be made concerning the 
Xmin values presented in Table in both approaches for 
Mb, the four values of Xmin a-re very close to each other. 
In particular, the differences between them are much less 
than the la interval (Ax'^ ~ 2.3). Therefore, as expected, 
the present data are not sufficient to discriminate between 
these four quintessence models even at the 1 a confidence 
level. 



5 QUINTESSENCE MODELS AND FUTURE 
SNEIA OBSERVATIONAL DATA 

The analysis of the present SNela data discussed in the 
previous section clearly showed that more observations are 
needed to clearly determine the nature of the dark energy. 
Such data could be provided by the proposed SNAP ( "Su- 
perNova/ Acceleration Probe") mission (see SNAP URL). 



SNAP, if accepted, will be a two-meter space telescope ded- 
icated to the discovery and observation of SNela with red- 
shifts in the range < z < 1.7. This instrument should 
perform photometry and spectroscopy of more than 2000 
SNela per year. The aim of this section is to investigate 
which values of the parameters (chosen by Mother Nature), 
future SNAP observations will be able to discriminate. The 
analysis will also reveal degeneracies between individual pa- 
rameters of the different models. 

We have first generated mock SNAP data by succes- 
sively assuming that model A, B, C or D was the true one, 
and under the ideal hypothesis that Mb was exactly known^. 
For each model, we have selected 100 couples of parameters 
restricted to physical regions surrounding contours given on 
Figure |S] (see Table This is why the wo parameter has 
been limited to [— 1;0]. The w\ parameter of model C has 
been limited to the interval [—1; 1], in order to avoid large 
absolute values of the EOS at z ^ 1. The parametrization 
adopted to generate data from model D means that we re- 
stricted ourselves to universes different from an Einstein- 
de Sitter one (M < 5) or to universes characterized by 

< 0.6 (/ > 0.4 or M > 2) (cf. Figure 

Then, for each value of the parameters, we have gener- 
ated three samples. For each of them we have computed the 
magnitude of 2100 SNela, with 2000 uniformly distributed 
in the redshift interval z G [0; 1.2] and 100 in the redshift in- 
terval z € [1.2; 1.7]. Gaussian errors on the magnitude have 
been added with a standard deviation of 0.15 mag, reflecting 
only the intrinsic spread of the SNela maximum luminosity 
(see SNAP URL). So for each of the four models, we are 
left with 3 X 100 sets of SNAP data, for which we know the 
parent model. Then, each sample has been independently 
fitted with the four models. The results are presented in 
Figures 171 1^ II II and 1131 where the parent is model A, B, C, 
D respectively. The parent model is in the right column and 
is successively compared to the other models. The squared 
symbols select the value of the parameters for which the 
two models are degenerated at the 2 a level, i.e. for which 
A Xmin < 6.17, whereas the filled dots represent the param- 
eters for which AXmin > 6.17 . Since we simulated noisy 
data, fitting the true model may not yield the exact value 
of the parameters. That is why the sets of models do not 
reproduce a regular grid in the right column of Figures |7| 
1^ nil and 1131 Moreover, different best fit positions can be 
associated to several simulations. 



5.1 Analysis of data generated with model A 

Figure|7|shows results obtained by analysing data generated 
from model A. It shows that it could be difficult to unam- 
biguously distinguish a model A with f2fc = from a model 
B, C or D at 2 (J level. Some misleading results can also be 
obtained from a model A with Q.k 7^ 0. This displays degen- 
eracies between the curvature and the EOS of dark energy : 
a FLRW model with spatial curvature {Q,k 7^ 0) and a cos- 
mological constant {w{z) = — 1) can be confused with a spa- 

^ Other works directly dealing with luminosity distance recon- 
struction instead of magnitude reconstruction assume that Mb is 
perfectly known (Huterer & Turner 1999, Saini et al. 2000, Chiba 
& Nakamura 2000, Gerke & Efstathiou 2002). 
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Model Values used for generate samples 



A 


!^m,i=0.1-(l-i)/5 




flj^j = -0.9 + 0.2 (j - 1) 


B 


nm,i = 0.05 + (j - i)/io 




woj = -0.05 - (j - 1) / 10 


C 


wo,i = -0.05 + (1 - i) /lO 




wij = -0.95 (2 j - 11)/ 9 


D 


/i = 0.4 + 1.6(i-l)/9 




Mj = 2 + (j - 1) / 3 



Table 3. Values of the parameters used for the SNAP data sim- 
ulations {i.j = 1,...,10). 



Wn 



n 

-0.2 - 
-0.4 
-0.6 - 
-0.8 - 
-1 



^ *f ** 

♦ 



-0.5 



0.5 



1.5 



r7_j].' ^ • ' 



1 

-0.2 ^ 
-0.4 4 
-0.6 - 
-0.8 



-1 ' 



-1 





-1 - 



um 



mm 




Figure 7. Positions of the Xmin ™ four models of the 300 
simulated data samples generated with model A. The filled circles 
(resp. empty squares) correspond to ^X^inx-A 6.17 (resp. 
^xL„,x-A<6.17)forX = B,C,D. 



tially flat model containing a dark fluid with a more general 
EOS (see also Figure |S| below) . On the other hand, points 
located on the edges of the considered parameter spaces cor- 
respond to model A values for which it is unlikely to obtain 
a good fit with a wrong model, i.e. for which reliable results 
can be expected. Let us also mention the following points ; 
— With model B, misleading results are obtained if the val- 
ues of (r2„i,r2A) are confined in the region inside the two 



0.5- 



W| 0- 



-0.5 



♦ ♦ 
♦♦ ♦ 



♦ ♦ • 

"I' 



♦ 



-1.5 



-0.5 





a. 



0.5 



1.5 



Figure 8. Correlations between individual parameters of model 
A and model C (when Ax^ ■ < 6.17). Upper plot - between 
Qrn — of the correct model A and wq of the wrong model C, 
Lower plot - between flj, of the correct model A and wi of the 
wrong model C. 



lines f2ft = and Qa ~ 0. 

— Ambiguous results between models A and C are expected 
if the true model is model A with {Qm, ^a) satisfaying these 
three conditions : Qrn+i^A < 1, Qm — ^A < 1 and Qm > 0.25. 

— In comparison with model D, most of the squares are lo- 
cated on the line flk ~ with Qrn > 0.5 in the parameter 
space of model A. This means that there is no degeneracy 
between the EOS of a PNGB fiuid and the curvature param- 
eter. In addition, more than half of the Xmin d positions are 
located in a region where M ~ 2 and / ~ 0.4, i.e. charac- 
terized by flm > and by a scalar field which is presently 
frozen or is just starting to evolve through the minimum of 
its potential (cf. Figure |5J. 

Finally, Figure |H| shows the correlations rising between 
individual parameters of model A and model C when both 
models are undistinguished at 2 a level (AXmin < 6.17). Lin- 
ear regressions through the points yield : 

wo I {nm. - Qa - 1) , (19) 



Wl 



3 



+ nA - 1 



(20) 



This demonstrates the degeneracy between the curvature of 
an FLRW model and the wi parameter of a model C. Note 
that by construction, model C with wo = —1 and wi = is 
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Figure 9. Positions of the Xmin ™ four models of the 300 
simulated data samples generated with model B. The filled circles 



(resp. empty squares) correspond to A Xmin 



> 6.17 (resp. 



Ax' 



< 6.17) for X = A,C,D. 



Figure 10. Upper plot - Correlation between the best fit values 
of Qrrt obtained with the correct model B and with the wrong 
model A, when Ay^ . < 6.17. Lower plot - Correlation between 

' ^min ^ 

the best fit values of wq obtained with the correct model B and 
with the wrong model C, when Ay^j^j < 6.17. 



identical to an FLRW model with Qm ~ 0.3 and JIa ~ 0.7. 
This result is basically reproduced by the relations (1191 and 



5.2 Analysis of data generated with model B 

Figure|^shows the expected degeneracies if the true model is 
a model B. The irregurar distribution of the best fit positions 
in the parameter space of model B denotes its stronger sen- 
sitivity to the statistical noise. Moreover, since we looked at 
the Xmin scanning parameter space from £7™ = to fi™ — 1, 
we find an accumulation of points in the upper left corner of 
the parameter space of model B, as expected from the form 
of the degeneracy curves shown in Figure |3] Note that the 
first row of Figure|^shows again the degeneracy between i}k 
and Wo, as already noted in Figure |7|: a fiat model B with 
Wo ^ —1 can be confused with an FLRW model with spatial 
curvature. 

Figure IIUI shows the possible error on the estimate of 
Qm and wq, when the correct model B is not used to fit the 
data. Although there is no significant difference in terms of 
Xmirn fitting the data with an FLRW model yields values for 
Sim larger than or equal to the exact value (upper plot), and 



fitting the data with a model C often give rise to value of 
Wo larger than or equal to the correct one (lower plot). 

5.3 Analysis of data generated with model C 

Figure 1111 shows the best fit positions in the parameter 
spaces of models A, B and D when model C has been used 
as a fiducial model. One can see that : 

— First, the Xmin, a &re all located in a region charaterized 
by ilm + < 1 and Qrn — ^a < 1- This confirms what 
has been obtained in Figure |7| (cf. second plot of model A 
on Figure |7|l. 

— Second, as expected, the Axmin between models B and 
C less than 6.17 are usually associated to the model C with 
wi ~ 0, but only if wo > —0.4. 

— Third, as previously noted, most of the Xmin positions 
in the parameter space of model D are located in a region 
where the scalar field is presently starting to evolve toward 
the minimum of its potential, i.e. Wx{z) leaves the value —1 
and approaches unity, with il^ > Qrn- This is why in most 
cases, the corresponding squares in model C are character- 
ized by wi ^ 0. 

Figure 1121 shows the correlation between the best fit 
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Figure 11. Positions of the X^in '"^ ^'^^ models of the 

300 simulated data samples generated with model C. The filled 
circles (resp. empty squares) correspond to ■^X^inX— C ^ 6.17 
(resp. A xLn,x-C < 6-17) for X = A,B,D. 



Figure 13. Positions of the . in the four models of the 

^min 

300 simulated data samples generated with model D. The filled 
circles (resp. empty squares) correspond to A x 
(resp. A x^i„ x-D < 6-17) for X = A,B,C 



Ln,x-D > 6.17 



1 
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Figure 12. Correlation between the best fit values of wi obtained 
with the correct model C and the best fit values of obtained 
with the wrong model B, when Ax^ ■ < 6.17. 



values of wi obtained with the true model (C) and the best 
fit values of flm obtained with a wrong model (B), when 
^Xmin < 6.17. If the true value of wi is positive (negative), 
then the value of fim found with the wrong model will be 
larger (smaller) than the true one (fim = 0.3). 



5.4 Analysis of data generated with model D 

Finally, Figure 1131 presents results obtained using samples 
generated with model D. It clearly appears that most of the 
Xmin positions corresponding to AXmin > 6.17 are located 
on the edges of the parameter spaces of the wrong models. 
So if the cosmological model chosen by Mother Nature is 
a PNGB model, it could be difficult to obtain a good fit 
with one of the three wrong models A, B or C. Moreover, 
the presence of some points outside the region of model D 
used for generating samples denotes a sensitivity to statis- 
tical noise. Note that most of the squares in the parameter 
space of model A of this figure appear like three bundles 
located at Jim ~ 1 — Qa ~ 0.3, 0.6 and 1. This is due to 
the strong dependence of Qm in model D with M when the 
scalar field is still frozen and acts like a cosmological con- 
stant (cf. horizontal lines for Qm in Figure |5J. The more 
precise is the grid used for the data processing in model D, 
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: -0.05 + (1 


-0/10 
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0.4 + 1.6( 
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f^A.j = 


= -0.9 


+ 0.2(i - 1) 
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- U - i)/io 
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- ll)/9 
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Compared with 


B 


C 


D 


A 


C 


D 


A 


B 


D 


A 


B 


C 


Pi^xLn < 11.8) 


19% 


31% 


14% 


98% 


93% 


82% 


95% 


68% 


40% 


35% 


23% 


23% 


^(^xLn < 6-17) 


15% 


27% 


9% 


95% 


93% 


72% 


79% 


59% 


26% 


29% 


20% 


21% 


Pi^xln, < -11.8) 


0% 


0% 


0% 


0% 


0% 


0% 


0% 


0% 


< 1% 


6% 


7% 


8% 



Table 4. Probabilities to obtain misleading results when the true model is model A, B, C and D at the 3 a level 
(P(Ax^i„ < 11.8)) and at the 2 a level (P(AxJ^i„ < 6.17)). P{Ax^i„ < -11.8) is the probability to exclude the right 
model at the 3 a level. All these values are associated to the limits on the parameter spaces explored for generating the 
data. 



the more those squares will appear uniformly distributed on 
the line flm + i^A = 1 in the parameter space of model A. 

An important remark has to be made at this point : 
Figures El 1111 and 1181 show that, whatever the true model 
is, a FLRW model with Sl„ = 0.3 and = 0.7 {w = ~1) 
cannot be rejected on the basis of the SNAP observations. 

5.5 Probability to select the correct model 

Table |1] shows the probabilities to obtain misleading results 
when the true model is model A, B, C and D at the 3 a 
level (P(Axmin < 11-8)) and at the 2a level (P(Axmin < 
6.17)). It also presents the probability to exclude the right 
model at the 3cr level (P(AXmin < -11.8)). They give an 
idea of the degree of degeneracy between the models, with 
respect to the computation of the luminosity distances. For 
example, it is difficult to identify a true model B among 
wrong models : the probabilities to confuse model B with 
model A, C and D are respectively 95 %, 93 % and 72 %, at 
2 a level. Model C cannot be distinguished from models A, 
B and D in 79%, 59% and 26% of the cases respectively 
{2a level). Note that the use of model C as fiducial model 
assumed an exactly known i},n. An uncertainty on Qrn will 
increase the percentages given above. Finally, if the fiducial 
model is A or D, we can expect to rule out the proposed 
wrong models with an a priori probability higher than 65 % 
at the 3 a level. However this does not preclude a priori a 
degeneracy between those models any other possible model. 

Samples simulated with model D can lead AXmin values 
less than —11.8 (~ 7 %), which means that the correct model 
can be excluded at the 3 a level and confirms the sensitivity 
of model D to the statistical noise. Table^Jalso shows that, 
whatever the true model is (B, C or D), the probability to 
reject model A at 3 a is low. 

The probabilities presented in Table |3] clearly depend 
on the regions of the different parameter spaces used for 
the data simulations. However, in order to obtain meaning- 
ful results, we chose the parameter ranges surrounding the 
physical part of the contours presented in Figure |S] 



6 CONCLUSIONS 

Making use of four quintessence models, of which three are 
mathematically simple and one is coming directly from par- 
ticle physics, we have first shown that the various cosmo- 
logical models may predict apparent magnitudes for objects 



at given redshifts which differ from each other by less than 
0.04 mag till z — 2. These magnitude differences are thus 
small compared to the intrinsic spread of the SNela maxi- 
mum luminosity (0.15 mag). This indicates that an unam- 
biguous discrimination between the cosmological models will 
be difficult to reach, with cosmological tests only based on 
luminosity distances. 

Second, we have fitted these models to the present su- 
pernovae data and found equally good fits with the four 
models. Finally, we have explored the discriminatory power 
of future SNela data to constrain the dark energy proper- 
ties. We have simulated a large number of SNAP data with 
each particular model and have re-analysed them with the 
four cosmological models. We have then compared the posi- 
tions and the values of the Xmin obtained using wrong mod- 
els with those obtained with the fiducial model. This led us 
to confirm the more skeptical conclusions already made in 
some previous works on the subject (Barger and Marfatia 
2001, Astier 2001, Maor et al. 2001, 2002) : some degen- 
eracies between the curvature, the matter density and the 
equation-of-state of dark energy will be difficult to break at 
the 3 a level and even at the 2 a and 1 a level. We also found 
that different estimates of a basic parameter like Qrn can be 
obtained, depending on the model used for the data pro- 
cessing. Moreover, whatever the true model is, the presently 
admitted FLRW model with Q„, = 0.3 and Qa = 0.7 will 
not be rejected on the basis of the future SNAP observations 
alone. 

In order to compare results expected from SNAP with 
the ones which could be obtained by other supernovae sur- 
veys, we have also generated data samples from models C 
and D for another redshift distribution : 2050 SNela between 
z = and z = 0.5 and 50 between z = 2 and z = 2.5 (with ai 
= 0.15 mag). This distribution reflects future observations 
to be made with the 4m Liquid Mirror Telescope (see ILMT 
URL) and with the Next Generation Space Telescope (see 
NGST URL). No significant difference has been found with 
the results presented here. 

Many recent papers discussed the strong constraints 
that could be expected from future SNela data on cosmo- 
logical parameters. However, the most important question 
is : how can we be sure that we consider the right model ? 
Indeed how useful are strong constraints on cosmological pa- 
rameters if they do not describe the model chosen by Mother 
Nature ? The degeneracy problem presented in this paper is 
expected to affect every cosmological test using luminosity 
distance. Therefore it will be very difficult for SNAP alone 
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to obtain any strong and safe constraints on the cosmolog- 
ical parameters. However, we expect that the combination 
of the CMB data with the future SNAP ones could help 
to break some of the luminosity distance degeneracies. This 
fact is already known in the framework of a single cosmo- 
logical model (see e.g. Huterer & Turner 2001, Frieman et 
al. 2002, Bean & Melchiorri 2002). 
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